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1.0 Introduction 

The work that follows is concerned with the application of design methodologies 
for feedback compensation associated with linear systems. In general, the intent is to 
provide a well behaved closed loop system in terms of stability and robustness (internal 
signals remain bounded with a certain amount of uncertainty) and simultaneously 
achieve an acceptable level of performance. The approach here has been to convert the 
closed loop system and control synthesis problem into the interpolation setting. The 
interpolation formulation then serves as our mathematical representation of the design 
process. Lifting techniques have been used to solve the corresponding interpolation and 
control synthesis problems. Several applications using this multiobjective design 
methodology have been included to show the effectiveness of these techniques. In 
particular, the mixed H 2 -H°° performance criteria with algorithm has been used on 
several examples including an F-18 HARV (High Angle of Attack Research Vehicle) 
for sensitivity performance. 


2.0 An Interpolation Viewpoint 

In this section, we will review some techniques that enable the control designer to 
view the closed loop control synthesis problem as an interpolation type problem. To 
see this, consider the problem of finding a compensator c to yield an acceptable level of 
closed loop performance and robustness (for a given open loop plant g) while 
maintaining internal stability. A typical block diagram for this problem is shown in the 
following Figure 2. 1 



Figure 2.1 Closed Loop System 

where ri and X 2 are exogenous signals (disturbances), y is the controlled output signal, 
and the error e and control input with noise u are internal signals. A compensator c is 
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called internally stabilizing if all of the elements of the following transfer matrix 


f 

e 


(l+cgr 1 

-g(l + eg) 1 

r > 

r i 

u 

L J 


c( 1 +Cg ) -1 

(1 +Cg) -1 

T2 

V ^ 


are in H°°. By H°°, we mean the set of all uniformly bounded analytic functions in the 
open unit disc. In engineering terms, H“ is simply the set of all stable transfer functions 
for discrete time systems. In particular, a rational function g is in H°° if and only if all 
of its’ poles are in the complex region | z | > 1 . 

Now for simplicity of presentation, assume that all the signals in Figure 2. 1 are 
scalar valued. Then in view of Figure 2.1, a typical controls synthesis problem may be 
to design an internally stabilizing compensator c to contain the effect of the 
disturbances rj and r 2 on the internal signals e and u or output signal y. For example, 
consider the signal to noise ratio (transfer function) from rj to e in (2.1) given by 


e _ 1 

r i 1 + eg ' 


( 2 . 2 ) 


Then setting a prespecified tolerance y, the control synthesis problem here is to find an 
internally stabilizing compensator c that norm bounds the transfer function 


11—11 = 11(1 +cg)-‘||<v (2.3) 

r l 

where the norm ||*|| is chosen by the designer. In the H°° literature, this is generally 
called bounding the sensitivity function, see for example [3,28,87]. The solution to this 
problem yields a stable closed loop system assuring that the signal to noise ratio is 
bounded by the inequality ||e|| < y ll r i II when the disturbance r 2 is neglected. Thus, 
setting y small decreases the effect of the external disturbance rj on the error signal e. 
As seen by the expression in (2.1), designing a compensator c to internally stabilize the 
system on all input-output transfer functions and satisfy a performance criteria as in 
(2.3) is nonlinear and awkward to work with. An alternate way of viewing the control 
synthesis problem is to first parameterize all stabilizing compensators in terms of one 
free parameter h in H°°, and then design h for particular closed loop performance 
objectives. This latter conversion technique enables us to view the control synthesis 
problem as an interpolation type problem. 

To show how the closed loop control system synthesis problem transcends into an 
interpolation problem, first assume for simplicity that g is rational. With this 
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assumption g admits a coprime factorization of the form g = n/d where n and d are 
rational functions in H°°, and there exists functions p and q in H°° that satisfy the 
following Bezout identity 

np + dq=l, {n,p,d,q} € H” . (2.4) 


Moreover, there exits several efficient algorithms to compute the H°° functions n, p, d, 
and q such that g = n/d and (2.4) is satisfied. Using the factorization g = n/d, it has been 
shown in [85] that the set of all internally stabilizing compensators for the closed loop 
system in Figure 2. 1 is given by 


p + dh 
q - nh 


he H~ 


(2.5) 


where h is an arbitrary function in H°°. It can be easily confirmed that substitution of c 
into (2.1) yields the internal transfer matrix 


e 

u 


(qd-ndh) (-nq+n 2 h) 
(dp-d 2 h) (qd-ndh) 


r i 

r 2 
v J 


( 2 . 6 ) 


Since n,d,p, and q are all in H°°, equation (2.6) shows that all the components in the 
transfer matrix in (2.1) are in H°° for arbitrary h in H°°. Therefore, all compensators c 
of the form (2.5) are internally stabilizing. One can also show that all internally 
stabilizing compensators are of the form c in (2.5). 

Now let us return to the control problem introduced in (2.3). That is, find an 
internally stabilizing compensator c such that the transfer function from rj to e given by 
(1 + gc) _1 is bounded in norm by y. By substituting the compensator c in (2.5) into the 
1-1 component of the transfer matrix in (2.1), we see that 


(1 + eg) 1 = dq - ndh = f - 0h 


(2.7) 


where f = dq and 0 = nd are in H°°. (Notice that using c in (2.5), every element of the 
transfer matrix in (2.1) has a similar form as in (2.7) for an arbitrary h in H“.) 
Therefore, the control synthesis problem in (2.3) and other control synthesis problems 
concerned with norm bounding the elements of the transfer matrix in (2.1) can be 
converted to the following Sarason type interpolation problem. 


Sarason Problem Given a particular f and 0 in H°° with a tolerance y > 0, find h in H 


that satisfies the following inequality 
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||f- 0h|| <y (he H°°) . 


( 2 . 8 ) 


This problem first appeared in the mathematical literature [68], with one of the first 
applications to control theory in [18]. From (2.8), one solves for a particular h in H°° 
that satisfies the bound ||f-0h||<y. Then from this resulting h, the corresponding 
internally stabilizing compensator c satisfying the inequality in (2.3) is given by 
c = (p+dh)/(q-nh). Of course, the solution to this interpolation problem depends upon 
the norm chosen. Some norms are inherently difficult to compute with. For practical 
and computational reasons, we will concentrate on solving interpolation and control 
problems using the H 2 and H°° norms. Recall that the H 2 norm (respectively, H°°) for a 
transfer function g analytic in | z | < 1 is defined by 


llglb = 


2k 


2n 


^ 1/2 


J lg(e ll )| 2 dt 

o 


and ||g||oo = sup |g(e ll )| . 
0 < t < 2 n 


(2.9) 


0 

To relate these norms to some of the results mentioned in the historical review, the H 
norm was basically used as the criteria for the development of the Weiner filter and the 
LQG optimization problem, and the H°° norm was associated with classical control 
theory (Bode [9], etc.) and more recent results (for example, see [18,19,20,87,88]) 
based upon frequency domain techniques. From a practical point of view, the H 
controllers yield nice closed loop performance properties and H°° controllers are robust 
to uncertainties in the system. In practice, one has to compromise performance for a 
certain amount of robustness in the synthesis problem. To accomplish this, the 
synthesis problem here consists of designing controllers for the closed loop system that 
yield an acceptable level of H 2 performance and inherits the H°° robustness. This leads 
to mixed H 2 -H°° type control problems. The typical approach in the existing literature 
is to solve these types of problems using Riccatti equation techniques. Here, we will use 
commutant lifting techniques to solve these problems. The lifting approach considered 
here enables us to consider the more general irrational interpolation problems and 
readily leads to efficient computational techniques in the rational case. 

To be explicit in defining control problems of concern here, lets consider again 
the problem of designing an internally stabilizing compensator c that norm bounds the 
sensitivity function (1 + gc) -1 . That is, search for an internally stabilizing compensator 
c that yields an acceptable level of performance (H 2 norm bound) and contains a 
prescribed amount of robustness (H°° norm bound) by satisfying the constraints 
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||(1 -t-gcr'lb <p 2 and ||(1 +gc)" 1 |U< poo (2.10) 

where p 2 and are prespecified H 2 and H°° norm bounds, respectively. From our 
previous analysis, we can convert the sensitivity problem into an equivalent Sarason 
type interpolation problem and search for an h in H°° satisfying the bounds 

||f- 0h|| 2 < p 2 and ||f-0h|U< Poo . (2.11) 

Thus, setting c = (p+dh)/(q-nh) from h in (2.11) (where g = n/d and the pair p and q 
satisfy the Bezout identity np+dq = 1) gives us a compensator c that satisfies the closed 
loop bounds in (2.10). 

In the sections that follow, we will assume that the reader is familiar with the 
Youla parameterization previously discussed as applied to a typical controls problem. 
We will first present some results obtained for the mixed H 2 and H°° control synthesis 
problem with examples including an application to the F-18 HARV model. Finally, in 
part II we will provide an algorithm for obtaining some mixed H 2 -H°° bounds for the 
continuous time control problem. 



Part I Discrete Time H 2 -H°° Compensator Design 
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1. Summary and Main Theorem 

Recall from the introduction that if we are given the plant g = n/d, (stable or not), 
we can parameterize all stabilizing controllers c by the expression c = (p+dh)/(q— nh) 
where n and d satisfy the Bezout identity pn + dq = 1 . This leads to the 
parameterization of all internally stable transfer functions given by various affine 
functions f- 0h where f and 0 are stable and h is an arbitrary function in H°°. With this 
in mind, we can view the sensitivity problem as the equivalent problem of norm 
bounding the error on a function of the form f - 0h where h is our design parameter. 
Now suppose we are given a specified function F in H°°(‘E ] , £ 2 ) (the multivariable 
setting) with the realization {A F , Bp, Cp} and {A©, B@, C©, D©} is a minimal 
realization for the inner function O(z). Let Q© F to be the solution to the following 
Lyapunov equation 

Q©f = AqQqfAf + C©C F . (1.1) 

Also let Wj be the controllability operator from \\(U) onto X generated by the pair 
{A©,C© } given by 

W, = C©, A©C©, A© 2 C©, ... (1.2) 

and let W2 be the controllability operator from l+(90 into X generated by the pair 
{A©,B} given by 

W 2 = B, A©B, A©B, ... (1.3) 

where B = Q©fB f . Define Pi and P 2 as the controllabiliy grammians for the pairs 
{A©,C©} and {A@,B}, respectively, where Pj and P 2 are the solutions to the following 
discrete time Lyapunov equations 

P, =W,Wl =A©P,A 0 + C©C© and P 2 = W 2 W 2 = A©P 2 A© + BB* . (1.4) 

Now define the H 2 and H°° optimization problems 

d 2 = inf { ||F — 0H|| 2 : H e H°° ) (1.5) 

d M = inf { ||F - 0H||oo : H e H“ } . (1.6) 

Then it is well known from the literature that d 2 can be computed from the relation 
d 2 =trace(B Pj'B) and di = ^ max (P7 1 P 2 ). Next choose 8>1 and set y=8d 00 . 
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Finally, let 

A p = (fPi - P 2 )A 0 (7 2 Pi - P 2 + BBV (1.7) 

B p = (7 2 Pi -P 2 )A 0 (y 2 P 1 -P 2 + BB*r'B 
C p = 7 2 C 0 ( T 2 P 1 -P 2 +BB*r' 

D p =7 2 C 0 (7 2 Pi -P 2 +Bb")- 1 B. 

be the finite dimensional operators. Then we have the following theorem. 

Theorem 1 

Let the finite dimensional operators Ap,B p ,C p , and D p be defined as in (1.7). 
Then A p has all of its eigenvalues inside the unit disc. Furthermore, letting G be the 
function defined by G(z) = D p + zC p (I - zA p ) _1 B p , there exists a function H in 
IT'X'E i , !f ) satisfying F - 0H = G. In addition, for this particular H, the following 
inequalities are satisfied 

5d 2 

||F - ©Hlloo < fid,, and ||F-0H|| 2 < (1.8) 

" dT 

1 "* T~ 

ndl 

where n is the dimension of ( ff. In particular, by choosing 8 2 = 2 — d 2 / nd 2 , F — 0H 
satisfies the following bounds 

||F - 0H|U < d_ -\/2 — ^4“ and ||F- 0H|| 2 < d 2 A [2 - (1.9) 

\ n&i y ndi 



The results presented in Theorem 1 are very useful in solving control synthesis 
problems in H°°, and this theorem guarantees that an H can be found to obtain the 
bounds given in (1.8). Furthermore, with a realization for G given in (1.7),H can be 
found by simply taking the adjoint of 0 to obtain 

H(z) = 0* ( l/Z)(F(z) - G(z)). (1.10) 

We will conclude this section by summarizing the algorithm for computing an H in the 
previous Theorem when F and 0 are given. 
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2. H 2 -H°° Discrete Time Algorithm 

In this section, we will summarize the algorithm to construct an H satisfying the 
norm bounds in (1.8) and apply the algorithm to several aircraft examples. 

(1). Compute a minimal realization {A©, B©, C©, D©} for 0. 

(1.. Compute a realization {Ap, Bp, Cp} for F. 

4c $ 

(2.. Solve the Lyapunov equation Q©p = A©Q©pAp + C©Cp 

(4) . Form the finite dimensional operator B = Q©pBp 

(5) . Compute the controllability grammians 

Pi = A©Pi A© + C©C© and 
P 2 = A©PzA© + BB . 

(6) . Compute d 2 = trace(B Pj 1 B) and d i = ^maxCPT 1 P 2 X 

(7) . Choose 8 > 1 and find the realization {A p , B p , C p , D p } for G. 

(8) . Solve for H(z) = ©‘(^(FCz) - G(z)) where 0*(^) = D© + B©(zl - A©) _1 C©. 

z z 


We will conclude this section with two examples; the first is a drone aircraft and the 
second is the NASA F-18 HARV (High Angle of Attack Research Vehicle) fighter 
aircraft. 
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Example 1 Drone Aircraft 


In this section we will show by example how the results in section 1 can be 
applied to a typical aircraft controls problem. We borrowed the open loop plant data 
from [19] for the lateral attitude dynamics of a drone aircraft. The open loop plant is 
given by the state space model 


x = Ax + B(u + w) , y = Cx 


( 2 . 1 ) 


where 


r 

-.0853 -0.0001 -0.9994 0.0414 

0 

0.1862 


f *\ 

0 

-46.86 -2.757 0.3896 

0 

-124.3 

128.6 


0 

-0.4248 -0.0622 -.0671 

0 

-8.792 

-20.46 

, B = 

0 

0 1 0.0523 

0 

0 

0 


0 

0 0 0 

0 

-20 

0 


20 

0 0 0 

0 

0 

-20 

> 


0 


C = 


1 0 0 0 0 0 


D = 0 


( 2 . 2 ) 


The states in this model represent, in order: sideslip angle, roll rate, yaw rate, roll angle, 
elevon surface deflection, rudder surface deflection, elevon servo command, and 
rudder servo command. For simplicity, we will consider the single input to single 
output map from the elevon servo command to sideslip angle, resulting in a non- 
minimal realization of order 6. The input output map of the open loop plant g is then 
given by the following transfer function 

176.09s 2 + 226.55s-10.47 (2 3) 

S ~ s 5 + 22.91s 4 + 58.21s 3 + 4.08s 2 + 74.21s + 2.67 ' 


This resulting transfer function is non minimum phase with two unstable complex 
poles. The following paragraph summarizes the results of our simulations using the 
H 2 -H°° algorithm in section 2 for the open loop plant g given in equation (2.3). 

To begin, lets consider again the H 2 -H°° sensitivity problem as presented in 
section 1 with the given open loop plant g in (2.3). Recall that we want to find an 
internally stabilizing compensator c that bounds the errors 
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11(1 +gc) 1 II 2 — 72 and 11(1 +gc) 'lU^Yoo (2.4) 

where 72 and 7,0 are bounds given as functions of the parameter 8 > 1 in equations (1.8). 
The first step in our design process was to discretize the given plant g in (2.3) using the 
Bilinear transformation in the Matlab subroutine. We then applied the algorithm from 
section 1 for the construction of the compensator c with four different values for 8. 
Using fast Fourier transform (FFT) techniques, the H 2 and H°° norms were computed 
for the resulting closed loop transfer functions for values 8= 1.0001, 1.05, and 1.1. 
Recall from section 1 that for each particular 8 > 1 chosen, the main Theorem 
guarantees that the closed loop transfer function obtains the given bounds in (1.8). In 
our example, the actual closed loop H 2 and H°° norms of the closed loop transfer 
v function were computed for each 8. 

For example, fo the three values of 8, the actual H 2 norm was computed and 
compared to the bounds as shown in the following Table 2.1. 


Table 2.1 (H-2 Norm) 


8 

inf||(l +gc) ’lb 

ll(l+gc) 1 II 2 

72 

1.0001 

0.4202 

1.0288 

1.0294 

1.05 

0.4202 

0.7484 

0.8507 

1.1 

0.4202 

0.6395 

0.7533 


Notice from the data that as 8 -> 1, the actual H 2 norm ||(1 + gc) _1 1| 2 increases because 
the theorem predicts H°° optimality as 8 — » 1 and hence no guarantee on the H 2 
performance. As 8 increases, for example, for 8= 1.1, Table 2.1 shows that the closed 
loop has better resulting H 2 performance. To show how the H~ performance varies 
with 8, we have listed values of ||(1 + gc) -1 || M as a function of the same three values of 
8 as shown in the following Table 2.2. 
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Table 2 (H-infinity Norm) 


8 

inf||(l + gc) MU 

llo+gcrMU 

Yoo 

1.0001 

1.0299 

1.0300 

1.0300 

1.05 

1.0299 

1.0777 

1.0814 

1.1 

1.0299 

1.1229 

1.1329 


The optimal and computed norms are listed to show that they actually fall within the 
predicted bounds given in equations (1.8). It is also interesting to note that the actual 
closed loop poles approach unity as 8 approaches unity (H°° optimal) although this is 
the design for least sensitivity to uncertainty. The spectral radius of the closed loop 
system and compensator order (complexity) is given in Table 2.3 below. 


Table 2.3 (Continuous Time) 


8 

Spectral 

Radius 

Controller 

Order 

No. Unstable 
Controller Poles 

1.0001 

0.9995 

11 

3 

1.05 

0.9976 

10 

3 

1.1 

0.9974 

10 

3 


As expected, the peak on the Bode plot of the closed loop transfer function (1 + gc) -1 
becomes smaller as 8 approaches 1 . In addition, the following plots in Figure 2.4 show 
that the infinity norm bound is successfully achieved. 





15 


Figure 2.4 


Magnitude Plots of inv(1+gc), delta=1 .0001,1.05,1.1 

i , , 1 1 r 


1 1 i i i i 1 1 J 

0 1 2 3 4 5 6 7 



Since the idea behind minimization of the sensitivity function (1 +gc) -1 using 
H°° as a performance measure is to .increase robustness, we have plotted the resulting 
magnitude plots of gc. By transforming the discrete time closed loop system to the 
continuous time via the bilinear transformation, the resulting Nyquist plots of gc in 
Figures 2.5 on the following page clearly varify the expected increase in robustness as 
5 — > 1. 






Figure 2.5 


Nyquist Plots of gc, Continuous, delta=1. 0001 ,1.05, 1,1 



The resulting gain margins were computed for each 8 chosen and the results are given 
in Table 2.6 below. 


Table 2.6 (Continuous Time) 


8 

Gain Margin 

1.0001 

3.13 

1.05 

2.70 

1.1 

1.16 


Finally, Figure 2.7 is a Bode plot of the resulting continuous time closed loop transfer 
function ( 1 + gc) -1 for the various values of 8 > 1 . 







Figure 2.7 


Cont. Time Bode Plots of inv(1 +gc), delta=1 .0001 ,1 .05,1 .1 



This plot varifies that the peak magnitude on the Bode plot is reduced with decreasing 8 
as anticipated. 

The results and example in this section have shown to be a successful 
methodology for multiple objective control system synthesis. In particular, if the 
designer is interested in compromising objectives with respect to the H 2 and H°° norms, 
the resulting algorithm is numerically pleasing and does not require a search for 
obtaining reasonable results. The next example is a direct application of these results to 
the F-18 HARV fighter aircraft model. 




Example 2 F-18 HARV Fighter Aircraft 


We have successfully demonstrated a methodology for the construction of 
compensators that achieve closed loop bounds on both the H 2 and H°° performance 
problems. The previous results extend to the multivariable case and to the continuous 
time problem. In this section, we will show how these results can be implemented on 
the F-18 model that NASA Dryden has generously provided. 

In the example that follows, we are given a 151 order state space model of the F- 
18 HARV (High Angle of Attack Research Vehicle) aircraft at flight conditions of 
Mach .9 with 10 inputs and 10 outputs. This model includes 102 states specifically for 
high fidelity actuator dynamics augmented with the modal data of the aircraft structure 
itself. The actuator dynamics eith order of complexity (dynamical order) are as follows; 


aileron 8 
stabilitator 1 1 
rudder 7 

leading edge flap 8 
trailing edge flap 9 
vanes 8 


Since there are two sets of actuators, we have 51x2 = 102 states for 12 actuators 
inclusive in the model. For the main plant model of the aircraft structure, we have the 
states <[>, theta, and \|/, for bank angle, pitch, and heading angle. For displaced degrees 
of freedom, we have u, v, w for forward, side, and vertical displacements. Finally, 
p, q, and r are the roll, pitch, and yaw degrees of freedom in that order. This model 
also includes 10 symmetric modes and 10 antisymmetric modes and their rates which 
make up 40 states. The interested reader can consult NASA TM 4493 for details 
concerning the actuator dynamics and the augmented state space model. In order to 
demonstrate the effectiveness of our compensator synthesis methodology, we have 
taken the input-output map from thrust vectoring in the pitch mode to the true angle of 
attack a. This model includes 8 states for the actuator. In the aircraft modal data we 
have kept the lowest four symetric and lowest four antisymetric modes. Hence, our 
plant model consists of 25 states from which we will apply our design methodology. 
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In the control synthesis problem that follows, our design goal is to achieve good 
closed loop mixed H 2 -H°° performance using the sensitivity function (i.e. (1+gc) -1 ) as 
our resulting closed loop input output map. Here we have applied the mixed H 2 -H°° 
control synthesis algorithm as listed in this section to this single input - single output 
plant. For various 5 > 1 we have constructed compensators to achieve closed loop H 2 
and H°° performance bounds as given in Theorem I. For values of 5 = 1.0001, 1.5, 3, 
and °o we have constructed a state space solution for the closed loop transfer function 
g c i and computed the resulting H 2 norms. The following Table 2.8 lists these values, 
and show in particular that as 5 — » °° we obtain the optimal H 2 solution as shown by the 
last row of numbers in Table 2.8. 


Table 2.8 H 2 Norm Data 


8 

inf ||(l+gc) 1 1| 2 

ll(i+gcr l ll 2 

li 

1.0001 

275.3202 

304.5828 

304.5921 

1.05 

297.8435 

301.4759 

3.0 

275.5405 

278.1628 

OO 

275.3202 

275.3203 


The bound y 2 as guaranteed from Theorem 1 is obviously satisfied for all three values 
of 8 > 0. We have also computed the resulting H°° norm of g c ] for four values of 8 and 
have listed these values in the following Table 2.9. 


Table 2.9 H°° Norm Data 


8 

inf ||(l+gc) 1 |U 

Il0+gc) 'lloo 

Joo 

1.0001 

304.5989 

304.6121 

304.6293 

1.05 

310.8247 

319.8288 

3.0 

371.3927 

913.7967 

OO 

381.0382 

OO 
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Notice that as 8 -> 1, the resulting bound y^ is the optimal H°° bound and is 
successfully achieved as the first row of Table 2.9 suggests. 

The resulting Bode plots for 5 = 1.0001, 1.05, and °° are shown in the following 
Figure 2.10. 



Figure 2.10 

Its clear from this plot that as 8 approaches 1, the resulting magnitude plot of 
lid + gc) 1 1| appears flat across the spectrum. In general, one can place weighting 
transfer functions in the optimization problem to shape these magnitude plots as desired 
for high frequency roll off, etc. 

In some problems, the time domain characteristics are more critical to the 
designer such as the impulse response of the resulting closed loop system. For pure 
interest, we have plotted the impulse response of the resulting closed loop sensitivity 
function for two values 8= 1.0001 and 8 = °° and superimposed them to show the 
difference in responses as shown in the following Figure 2. 1 1. 
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Impulse Response for del = [1.0001 inf] 



Figure 2. 1 1 


As one would expect, the H 2 optimal solution (8 = °°) has a better damped response, 
but in general the H 2 is much less robust to uncertainty. 

In general, one dissadvantage of using interpolation type results for control 
design methodologies is that the resulting controller is in general high order. This can 
cause implementation problems due to the complexity of the controller structure. In this 
example, the resulting compensator was typically of order greater than 70-80 states, 
which makes it difficult for implementation. We had success in computing these 
compensators due to the numerical efficiency using fast Fourier transform techniques. 
For multivariable applications, one can use the discrete Riccatti equation for obtaining 
the compensator c via the Youla parameterization. 

In conclusion, this example was a more realistic problem that one encounters in 
industry. We have shown for the F-18 HARV single input single output system from 
pitch vectoring to angle of attack that the algorithm for mixed H 2 and H°° bounds has 
successfully achieved these bounds anticipated by Theorem 1 . 
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Part II Continuous Time H 2 -H°° Compensator Design 

/ 
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1. Summary and Main Theorem 

In the previous section, solutions were presented that solved variations the 
Sarason optimization problem. Performance was measured in terms of the H 2 and H°° 
norms. Thus far, everything has been presented in terms of the H and H°° spaces on 
the unit disc. In this section, we will solve these same types of problems on the spaces 
associated with the H 2 and H°° problems on the half planes. It can be shown that an 
appropriate unitary bilinear transformation from the unit disc to the half plane indeed 
allows us to consider similar H 2 - H°° type interpolation problems as done in section I. 
In particular, we will solve the continuous time Sarason problem and present an 
example to show how the algorithm we have developed can be applied to a typical 
controls problem. 


2. H 2 and I I 00 Optimization Problems 

In this section, we have made use of the commutant lifting theorem to extend and 
tighten the Kaftal-Larson-Wise bounds [47] to the right half plane setting. Throughout, 
L 2 (£) is the space of all square integrable Lebesgue measurable functions f(iw) for 
-oo<w<°o with values in the Hilbert space £, under the inner product 
(f,g) = l/(2rt)J(f,g)dw. The space H 2 (£) is the Hardy space of all analytic functions in 
the right half s-plane (Re(s) > 0) with values in £ whose norm 

oo 

I |f| 1 2 = supy- J ||f(a+iw)|| 2 dw (1.1) 

^ — co 

O>0 

is finite. By passing limits to the iw axis a.e., H 2 (£) is viewed as a subspace of L 2 (£). 
Throughout, K 2 (£) is the orthogonal complement of H 2 (£) in L 2 (£). If £ is the 
Hilbert space generated by the set of all Hilbert-Schmitt operators from U into y, then 
L 2 (£) is denoted by L 2 (U,90 and H 2 (£) is denoted by H 2 (?/,(X)- As expected, 
L 90 is the Banach space of all uniformly essentially bounded Lebesgue 
measurable functions on the interval (-<x>,oo) whose values are linear operators from U 
into y. The space H°°(£i ,90 is the subspace of ,90 corresponding to the set of 
all uniformly bounded analytic functions in Re(s) > 0 whose values are linear operators 
from U into y. 



24 


To expedite the main results of this section, we will summarize the proof of the 
optimal H 2 solution and state the results of the optimal H°° solution. To begin, let a be 
a scalar in C 1 with Re(a) > 0 and let <J) be the function in H 2 denoted by 
<|> = Re(a)/(s+a). Also let ||-|| 2 denote the L 2 (£Z,90 norm and assume U is finite 
dimensional. Let F be a specified function in L °°(W,90, and 0 an inner function in 
H“(£,90, that is, 0(iw) is a.e. an isometry. Now consider the weighted H 2 
optimization problem, 

d 2 ,o = inf{||(F-©H)(l»|| 2 :He H°°(?/,£)}. (1.2) 

For the solution to this problem, let P#' denote the orthogonal projection onto the 
space Of* = L 2 (90 ©■ 0H 2 (£). Then an application of the projection theorem yields 
d 2,a = llP^'F^lb- 

Now consider the H” Sarason optimization problem; 

d oo =inf{||F-0H|U:He H~(«,£)} (1.3) 

where F is a specified function in L”(£/,90 and 0 is an inner function in H°°(!E,90 
(all solutions to a more general 4-block problem are given in [36]). By using standard 
commutant lifting techniques [4], it follows that the optimal error is given by d^ = ||r|| 
where T is the operator from H 2 (U) into !H' defined by T = P^'Mp |H 2 ( ( U). (Here, 
Mg denotes the multiplication operator from L 2 (£i) into L 2 (90 defined by M G u = Gu 
for u in L 2 ( U ) where G is in L°°(£/ ,90). In fact, this result will be a by product of our 
analysis. The optimization problems in (1.2) and (1.3) occur naturally in engineering 
applications and control theory , see [4,8,28,39]. 

The above nomenclature and identifications enable us to solve a right half plane 
version with respect to the main theorem of section I in concern with H 2 -H“ 
optimization problems. Here, we have assumed that the typical controls problem has 
been converted to a Sarason type interpolation problem. Now we can state the main 
results of this section. 

Theorem 2 Let F be an element ofU°(U ,90. S > 1, and n the dimension of U . Also 
let 0 be an inner function in H°°(£,90. Then there exists an H a in H ( c ll, r E) 
satisfying ||F - ©H^l^ < Sd^ and 

( o l -l/2 

||(F - 0H a )<t>|| 2 < bd 2 a 8 2 -l+2d 2>ct /(nRe(a)di) 


(1-4) 
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In particular, by choosing 5 2 = 2-2d 2>a /(nRe(a)di), there exists an H a satisfying the 
following bounds 


||(F - 0H a )(J)|| 2 < d 2 , a \2 - 2d 2 a /(nRe(a)di) 
||(F - 0H a )|U < d J 2 - 2d 2 a /(nRe(a)di) 


1/2 

1/2 


(1.5) 

( 1 . 6 ) 


Notice in the previous Theorem 2 that our results depend upon the a chosen. 
This a can be used as a weight in typical H 2 and H°° controls problems. Our proof of 
Theorem 2 is based upon the central solution to the commutant lifting theorem 
[26,27,73,74]. The following paragraph discusses implementation of the results in 
Theorem 2. 

To show the practical aspects of Theorem 2, we will present some explicit 
formulas that enable us to construct an H satisfying the inequalities in (2.4). Finally, 
we will present some examples of of these results. 

We say that F is a rational function in K°°(U ,90 if F(s) is a proper rational 
function with all of its poles in the open right half complex plane. Let A be an operator 
on a finite dimensional vector space X, B the operator from U into X, and C the 
operator from X into y. Then {A,B,C} is a minimal anticausel realization for F if the 
pair {A,B} is controllable, {C,A} is observable, and 

F(s) = C(-sI-A) -1 B, (2.7) 

(see [48,66] for further details on realization theory). Since F is in ,y ), all of the 

eigenvalues of A are in the open left half plane. Since 0 is two sided inner, we can 
assume without loss of generality that 0 = 1, hence the operator T previously defined is 
now the Hankel operator mapping H 2 ( Z1 ) into 9f = K 2 (90 = L 2 (90 & H 2 (90 defined 
by r = P#'M F | H 2 (Z1), where Py/' is the orthogonal projection onto the space K 2 (90- 
It is well known (see [26,28,35,39,88]) that this Hankel operator T admits a 
decomposition of the form T = W 0 W* , where W 0 is the operator from X into K 2 (90 
defined by W 0 = C(-sI-A) -1 and W c is the operator from X to H 2 (7i) defined by 
W c = B*(sl-A*) _1 . The controllability grammian P and the observability grammian Q 
are defined by P = W*W C and Q = W*W 0 . Because {A,B,C} is a minimal realization, 
both P and Q are strictly positive definite operators and are the unique solutions to the 
following Lyapunov equations 

AP+PA*+BB* = 0 and A*Q+QA+C*C = 0 . 


(2.8) 
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The minimal realization {A,B,C} along with its controllability and observability 
grammians P and Q, will play a key role in presenting state space solutions for the 
function F-H a . Finally, recall that ||F|| 2 = d 2 equals the largest eigenvalue of QP, 
(see [4,26,28,35,39,88]). 

Using the previous definitions, the H 2 optimization problem is solved trivially by 
decomposing FO into orthogonal components. Notice from the following expansion, 

FO = C(-sI-A)" 1 B(Re(a)/(s+a)) (2.9) 

r 

= Re(a) C(al-A) -1 (s+a) _1 B+C(-s-A) _1 (al-A) -1 B , 

that the first and second terms of (2.9) are contained in H 2 (£Z,90 and K 2 ( < U,Q r ), 
respectively. So if 

H = C(al-A) -1 B , (2.10) 

then H is the solution to the following optimization problem 

d 2 ,a = inf{||(F-H)<t)||:H€ H°°} (2.11) 

= { ||(F - H)<|>|| : H = C(al - A) -1 B } . 

Moreover, the resulting optimal H 2 error is 

d 2 itt = ||Re(a)W 0 (aI-A) _1 B|| 2 = trace|^B*(aI-A*r 1 Q(aI-Ar 1 B (2.12) 

Since F is rational, a realization for G r = F - ©H a can be constructed as shown in 
the following Theorem 3. 

Theorem 3 Let {A on X,B,C} be a minimal anticausel realization for the function F 
in K~(£Z,90. P and Q be the Lyapunov solutions defined in (2.8) for the pairs 
{A,B} and {A,C}, respectively. Set 8 > 1 and y=5d 00 where di = ^max(QP)- N r 
in K°°( Zl , 90 and D r in H°°( Zl , U ) be the rational functions given by 

N r = y 2 W 0 (y 2 I-PQ) _1 (al-A) -1 B and (2.13) 

D r = — + WcQ^I-PQ) -1 (al-A) -1 B . (2. 14) 

s+a 

Then G r = N r D7' is a function in L °°(£/,90 satisfying the bounds in Theorem 2 and 
admits a decomposition of the form G r = F-H a for some H a in H°°(?i ,f). 
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Notice if we let a— » <» along the real axis, then (2.13) and (2.14) show that 
G r = N r D7' approaches G = ND -1 where 

N = y 2 W 0 (y 2 1-PQ) - 1 B and D =I + W c Q( 7 2 I-PQr 1 B . (2.15) 

Moreover, the function G = ND 1 = (F-H) for H in K°{ ( ll, < y) solves the Sarason 
problem with the bounds 

8||p^f|| 2 

||G|U < &L and ||G|| 2 < (2.16) 

Vs 2 — i 

for 5 > 1 . We conclude this section with a remark concerning the factorization of H 
satisfying G = F - H when 0 = 1. 


Remark Assuming a is positive real, the identity 

(sI-A*)<t> = (sI-A*)a/(s+a) = a[l-(s+a) _1 (aI+A*) 


(2.17) 


along with some tedious manipulations can be used to show that H a = F - G r has a 
realization H a = D a + C a (sl-A a ) -1 B a where 

A a = -{(y 2 Q _1 -P)+(aI-A) _1 BB* } _1 {(I-A/a) _l BB*-(y 2 Q _1 -P)A* } (2.18) 

B a = {(y 2 Q _1 -P)+(aI-A) _l BB* } _1 (aI-A) _1 B 

_ _ ^2 rn -i J{(7 2 Q' 1 -P>KaI-Ar 1 BB* I" 1 {(I-A/a)- , BB*-(7 2 Q- 1 -P)A* } 
a ~ r Q [ + (I-y- 2 QP)A*- Y - 2 aQP 

D a =7 2 CQ“ 1 {(VQ 1- P) + (otI-A) _1 BB* } -1 (al-A) -1 B . 


Furthermore, setting 

A h =A*-(t 2 I-QP)- 1 QBB\ B h = ("^I-QP^QB and C H =CP, (2.19) 

it can be easily shown that as a — > H a — » H = C H (sI-A H r , B H uniformly over the 

complex right half plane. As a matter of fact, this is the solution corresponding to the 
unweighted H 2 -H°° problem due to the fact that as a — > ot/(s + a) — > 1 . Notice if F 
is purely anticausel, the H satisfying the optimal H 2 unweighted problem (y,oc — > °°) is 
the trivial solution H = 0. The following section shows how these formulas are useful 
for rational functions with an application to control system synthesis. 


2. H 2 -H°° Discrete Time Algorithm 
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In this section, we will summarize the algorithm to construct an H a satisfying the 
norm bounds in (2.4) and apply the algorithm to an example. 


(1) . Compute a minimal realization {Ap, Bp, Cp} for F. 

(2) . Compute the Lyapunov equations for P and Q in (2.8) 

(3) . Compute d a and d^ 

(4) . Set 5 > 1 and a > 0 

(5) . Compute {A a , B a , C a , D a } for H a 

We will now apply the previous algorithm to a typical controls problem 
associated with disturbance attenuation. Suppose we are given an open loop plant g and 
consider the design of a filter c that achieves a preset performance criteria on the 
resulting closed loop system. The following block diagram shows pictorially the plant 
g, filter c, and exogenous disturbance w and controlled output signal z. 


W 



Figure 6.3 Closed Loop System 


From the diagram, the input-output map from w to z is given by 


— (s) = gel 

w 


g 

1 +gc 


(6.3) 


For simplicity, we will assume g = n/d is rational and stable with an inner-outer 
factorization given by g = gjg 0 . With open loop stability in mind, a parameterization of 
all stabilizing controllers is given by c = h/(l - gh) for arbitrary h in H°°. Using this 
expression for c and substituting into the closed loop expression in (6.3) we obtain the 
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following norm expressions for the closed loop system g c |. 


llgcilk- = ll --- ;, f - — lk~ = ||g - g 2 h|k~ = l|g,*go - goh|| 2 ,oo • (6.4) 

1 + g(h/(l -gh)) 

Thus, setting g 2 h = P+(g*g 0 ) + h\ we can equivalently consider the control problem in 
the Nehari setting with a search for h in H°° and norms 

llgci Ik- = H p +(gi* go) + P-(g* go) " goh|k~ = ll p -(g* go) “ h'lkoo • (6.5) 


Now consider the following stable open loop plant g with inner-outer 
factorization g = g,g 0 given by 


g(s) = 


Wn(s 2 -2C n w n s+w^) 

(s+l)(s+2) 


s 2 -2C n w n s+w^ Wn(s 2 +2^ n W n S+W 2 ) 1 _ 
s 2 +2£ n w n s+w 2 ^ (s+l)(s+2) 


(6.6) 


Using various values for 8 and a with = .005 and w n = 10, we have computed the 
rational function P_(g* g 0 ) - h and simulated the Bode plots. 

The magnitudes in decibels for P-(gi*g 0 ) _ h (which has the identical Bode plot as 
g c i) with values of 8 = 1.001, 8 = V 2 , and 8 = °° with various a > 0 have been plotted 
and shown in the following figures. 







31 


The following Table 6.5 gives the numerical values with a= 10 from the resulting 
computations. 


J 


Table 6.5 Norm Data 


8 

a 

deo 

IlGcilL 

INI— bound 

d2,a 

IlGcA 

IN| 2 bound 

1.001 



1.9599 

1.9599 


3.9819 

4.1764 

<2 

10 

1.9579 

2.4887 

2.7689 

0.6159 

0.7088 

0.8626 

oo 



2.8677 

OO 


0.6159 

0.6159 


Notice that for each 8 chosen the resulting norms are within the bounds given from the 
theoretical results. The following paragraph discusses the computations required to 
obtain the compensator c. 

Recall that the compensator c is given by c = h/(l-gh). Hence, using the 
expression g 2 h = P+(g* g 0 )+h in c we obtain 

P+(g*go) + h 

c = — ; — ■ (6.7) 

go “ g(P+(gi go) + h ) 


In particular, we computed the compensator using the formula in (6.7) and the resulting 
closed loop transfer function for the cases that 8 = V 2 and a = 0.1, 10 and 1000 in the 
previous example. For instance, setting 8 = ^2 and a = 10, the resulting compensator 
was computed from the formula in (6.7) and given by 

-124.8s 2 -50.76s-12357.22 ON 

c = . (6.8) 

s 2 -44.83s-548.25 


It is of particular interest that the resulting closed loop g c ] computed from the actual c 
in (6.7) for 8 = V2 and a = 0.1, 10, and 1000 has the expected Bode plot as shown in 
Figure 6.6. 
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Bode Plot, gel, delta = sqrt(2) 



Figure 6.6 Magnitude Plot 

For applications and actual implementation purposes, it is usually desirable to 
have low order compensators for minimizing computation time. We have listed the 
dynamical orders of the compensator c and associated closed loop system g c ] and the 
resulting number of unstable compensator poles in the following Table 6.7. 


Table 6.7 Compensator Data 


a(8=\ / 2) 

compensator 

order 

closed loop 
plant order 

no. of unstable 
compensator poles 

0.1 

3 

5 

1 

10 

2 

4 

1 

1000 

2 

4 

0 


Notice from Table 2.5 that for low values of a, the resulting compensator has one 
unstable pole, and as a — > «> the compensator is stable. 

Conclusions 

The previous theorems and examples show that the multiobjective H 2 -H°° 
control problem can be viewed from a Sarason (model matching) point of view. Here, 
we have used a lifting approach to obtain our results once in the interpolation setting. 
State space algorithms have been prescribed for the construction of these solutions and 
resulting compensator. We have also shown that our solution satisfies some interesting 
H 2 -H°° bounds. Future work would involve generalizing these results to 2-block and 
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4-block settings (we have implicitly solved the 1 -block problem). Finally, we hve 
successfully applied these algorithms to the F-18 HARV fighter aircraft in a design 
example. 
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